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Abstract 

The lifetimes of the B° and B~ mesons are extracted from 29.1 fb _1 of data collected with the 
Belle detector at the KEK S-factory. A fit to the decay length differences of neutral and charged B 
meson pairs, measured in events where one of the B mesons is fully reconstructed in several hadronic 
modes, yields t w = 1.554±0.030(stat)±0.019(syst) ps, r B - = 1.695±0.026(stat)±0.015(syst) ps, 
and t b -/t w = 1.091 ± 0.023(stat) ± 0.014(syst). 

PACS numbers: 13.25.Hw, 12.39.Hg 
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Lifetime measurements of B and B~ provide key input to the determination of the 
Kobayashi-Maskawa matrix element |V^,| jlj. Moreover, the ratio of the lifetimes is sensitive 
to effects beyond the spectator model. In the framework of heavy quark expansion, the B 
lifetime ratio is predicted to be equal to one, up to small corrections proportional to 1/ mf [|J , 
where mj is the mass of the b quark. A recent theoretical study predicts tb-/t b0 = 1.07 ± 
0.03 [f|, which agrees with the world average of measurements, T B -/r-g = 1.073 ±0.027 [|J. 
The BABAR collaboration recently reported r B -/r w = 1.082 ± 0.026 ± 0.012 §. Accurate 
determinations of the B° and B~ lifetimes also provide essential input to analyses of CP 
violation and neutral B mixing. 

The analysis described here is based on a 29.1 fb~ data sample, which contains 31.3 x 10 6 
BB pairs, collected with the Belle detector || at the asymmetric-energy KEKB collider 
(8.0 GeV electrons against 3.5 GeV positrons) 0. Electron-positron annihilations produce 
T(4S') mesons moving along the electron beamline (z axis) with (ftrfr = 0.425, and decaying 
into B°B° or B + B~ . Since the B mesons are nearly at rest in the T(4S') center of mass 
system (cms), B lifetimes can be determined from the separation in z between the two B 
decay vertices. The average separation is ctb(0j)t ~ 200 fim, where r B is the B meson 
lifetime. In this analysis, one of the B decay vertices is determined from a fully reconstructed 
B meson, while the other is determined inclusively using the rest of the tracks in the event. 

The Belle detector consists of a three-layer silicon vertex detector (SVD), a 50-layer cen- 
tral drift chamber (CDC), an array of 1188 aerogel Cherenkov counters (ACC), 128 time-of- 
flight (TOF) scintillation counters, an electromagnetic calorimeter containing 8736 CsI(Tl) 
crystals (ECL), and 14 layers of 4.7-cm-thick iron plates interleaved with a system of resistive 
plate counters (KLM). All subdetectors except the KLM are located inside a 3.4 m diameter 
superconducting solenoid which provides a 1.5 T magnetic field. The impact parameter 
resolutions for charged tracks are measured to be a 2 y = (19) 2 + (50/(p/9sin 3 ' 2 6)) 2 fim 2 
in the plane perpendicular to the z axis and a 2 = (36) 2 + (42/ (p/? sin 5 / 2 6)f fim along 
the z axis, where f3 = pc/E, p and E are the momentum (GeV/c) and energy (GeV), 
and 9 is the polar angle from the z axis. The transverse momentum resolution is 
(a p J Pt ) 2 = (0.0030//5) 2 + (0.0019p t ) 2 , where p t is in GeV/c. _ 

B° and B~ mesons are fully reconstructed in the following decay modes: B° — > D + tt~, 
D*+n-, D*+p-, J/tpK° s , J/^jK* , B~ -> D n~, and J/if)K~ @. 

Charged pion and kaon candidates are required to satisfy selection criteria based on 
particle-identification likelihood functions derived from specific ionization (dE/dx) in the 
CDC, time-of-flight, and the response of the ACC. Electrons are identified using a combi- 
nation of dE/dx, the response of the ACC, and the position, shape, and total energy (i.e., 
E/p) of their associated ECL showers. Muon candidates must penetrate the iron plates of 
the KLM in a manner consistent with the muon hypothesis. 

The reconstruction and selection criteria for Kg — > 7r + 7r~ and J/ip — > £ + £~ (£ = e,/i) are 
described elsewhere || . K *° candidates are reconstructed in the K *° — > K~tt + channel and 
required to have an invariant mass within 75 MeV/c 2 of the average K*° mass. 

Photon candidates are defined as isolated ECL clusters of more than 20 MeV that are 
not matched to any charged track. tt° candidates are reconstructed from pairs of pho- 
ton candidates with invariant masses between 124 and 146 MeV/c 2 . A mass-constrained 
fit is performed to improve the 7r° momentum resolution. A minimum 7r° momentum of 
0.2 GeV/c is required. p~ candidates are selected as ix~ix° pairs having invariant masses 
within 150 MeV/c 2 of the average p~ mass. 

Neutral and charged D candidates are reconstructed in the following channels: D° — > 
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K~7i + , K~7t + 7t°, K~7t + 7t + 7t~ , and D + — > K _ 7r + 7r + . D* + — > D°7r + candidates are formed 
by combining a D candidate with a slow and positively charged track, for which no particle 
identification is required. 

To reduce continuum background, a selection based on the ratio of the second to zeroth 
Fox- Wolfram moments JTIJ and the angle between the thrust axes of the reconstructed and 
associated B mesons is applied mode by mode. 

The decay vertices of the two B mesons in each event are fitted using tracks that are 
associated with at least one SVD hit in the r-<p plane and at least two SVD hits in the r-z 
plane under the constraint that they are consistent with the interaction point (IP) profile, 
smeared in the r-<p plane by 21 jum, to account for the transverse B decay length. The IP 
profile is described as a three-dimensional Gaussian, the parameters of which are determined 
in each run (every 60,000 events in case of the mean position) using hadronic events. The 
size of the IP region is typically o x ~ 100 jum, a y ~ 5 /im, and o z ~ 3 mm, where x and y 
denote horizontal and vertical directions, respectively. 

In the case of a fully reconstructed B — > D^*'X decay, the B decay point is obtained 
from the vertex position and momentum vector of the reconstructed D meson and a track 
other than the slow tt + candidate from D* + decay. For a fully reconstructed B — > J/ipX 
decay, the B vertex is determined using lepton tracks from the J/if>. The decay vertex of 
the associated B meson is determined from tracks not assigned to the fully reconstructed B 
meson; however, poorly reconstructed tracks (with a longitudinal position error in excess of 
500 /im) as well as tracks likely to come from K s decays (forming the K° s mass with another 
track, or more than 500 fim away from the fully reconstructed B vertex in the r-<f> plane) 
are not used. 

The quality of the fit is assessed only in the z direction (because of the tight IP constraint 

in the transverse plane), using the following variable £ = (l/2n) Y% [( z l&ei ~ before V^before] 2 ' 
where n is the number of tracks used in the fit, -2 before and z after are the z positions of each 
track (at the closest approach to the origin) before and after the vertex fit, respectively, and 
^before i s the error of z before . A Monte Carlo (MC) study shows that £ does not depend on 
the B decay length. We require £ < 100 to eliminate poorly reconstructed vertices. About 
3% of the fully reconstructed and 1% of the associated B decay vertices are rejected. 

The proper-time difference between the fully-reconstructed and the associated B decays, 
At = t rec — t asc , is calculated as At = (z rec — ^ asc )/[c(/37)x], where z Tec and z asc are the 
z coordinates of the fully-reconstructed and associated B decay vertices, respectively. We 
reject a small fraction (~0.2%) of the events by requiring |At| < 70 ps (~45r^). The final 
event selection is based on requirements on the energy difference AE = — E™^ and 
the beam-energy constrained mass M bc = \/(E^ m ) 2 — (p c ™) 2 , where E^ m , E^, and p c ™ 
are the beam energy, the energy, and the momentum of the fully reconstructed B candidate 
in the cms, respectively. If more than one fully reconstructed B candidate is found in the 
same event, the one with the best x 2 f° r and M bc (and the invariant mass of the D 
candidate for the B — > D^*'X case) is chosen. For each channel, a square signal region 
is defined in the Ai?-M bc plane, corresponding to ±3cr windows centered on the expected 
means for AE and M bc . A standard deviation is about 3 MeV/c 2 for M bc and 10-30 MeV 
for AE depending on the decay mode. After all vertexing and selection requirements, we 
find 7863 B° and 12047 B~ events in the AE-M hc signal region. Figure shows the M bc 
distributions for all the B° and B~ candidates found in the AE signal region. 

We extract lifetime information using an unbinned maximum likelihood fit to the observed 
At distributions. We maximize the likelihood function L = Y\ i P(Ati), where -P(At) is the 
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FIG. 1: The beam-energy constrained mass distributions of the fully reconstructed B candidates 
(neutral on top, charged on bottom). The dashed curves show the background contributions. 



probability density function (PDF) for the proper-time difference At for each event and the 
product is taken over all selected events. The function P(At), expressed as 

P(At) = (1 - / ol ) [/si gJ Psi g (At) + (1 - / sig )P bkg (At)] 

+/olPol(At), (1) 

contains contributions from the signal and the background (P S i g and Pbkg), where P S i g is 
described as the convolution of a true PDF (P S i g ) with a resolution function (P S i g ) and Pbk g 
is expressed in a similar way: 

/+oo 
d(At')P J (At / )Pi(At- At'), (2) 
-co 

where j = sig, bkg. To account for a small number of events that give large At in both 
the signal and background (outlier components), we add a third component: P Q i(At) = 
G(At;cr i), where G(t;a) = exp(— t 2 / (2a 2 )) / (\/2ttct) is the Gaussian function. The global 
fraction (/ Q i) and the width (a \) are left as free parameters in the fit. 

The signal purity (f S i g ) is determined on an event-by-event basis as a function of AE and 
Mb c . The two-dimensional distribution of these variables including the sideband region is 
fitted with the sum of a Gaussian signal function P S i g (AP, Mt, c ) and a background function 
Pbk g (AP, Mb c ), represented as an ARGUS background function JTTJ in M^ c and a first-order 
polynomial in AE. The fraction of signal is then obtained as / S i g (AP, Mb c ) = P S ig/[Psig + 
Pbkg]- The true PDF for the signal is given by 

P sig (At) = -^expf-^V (3) 



2r B V t b 
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where t b is, depending on the reconstructed mode in the event, either the B° or the B~ 
lifetime. The resolution function of the signal is constructed as the convolution of four 
different contributions: the detector resolutions on z Tec and z asc (P rec and P asc ), an additional 
smearing on z asc due to the inclusion of tracks which do not originate from the associated B 
vertex (R np ), mostly due to charm and K~s decays, and the kinematic approximation that 
the B mesons are at rest in the cms (Rk). 

For a vertex obtained from multiple tracks, R mc and P asc are parameterized as a single 
Gaussian function, G (i; (s° + s}£j)cj-J (j = rec, asc) where <x,- is the proper-time error esti- 
mated vertex-by- vertex from the Zj error, and (s° + sj£j) is a scale factor based on the vertex 
fit quality £ defined above. The parameters for the scale factor, s°' , are determined from 
the lifetime fit. The typical value of the scale factor is 1.32 for R vec and 1.20 for P asc - For 
a vertex obtained from a single track fl2|], we use (1 — / ta ii)G (t; s m ain0j) + ft&nG (t; Stance), 
where s ma i n and s ta u are global scaling factors. The shape of P np is determined from MC 
data sample, separately for B° and B~ events. R^ is calculated analytically as a function 
of E B ms and cos 8 c B ms from the kinematics of the T(4S) two-body decay, where ^^ ns is the 
polar angle of the reconstructed B in the cms. The resulting At resolution for the signal is 
~1.56 ps (rms). 

The background PDF, Pbk g (At), is modeled as a sum of exponential and prompt 
components convoluted with R hkg (At) = (1 - f^f)G (At; s^ k J nV / ov^ T +o a ^?j + 

f^fG (At; sSfVoVec 2 + Oasc 2 )- Different values are used for s bkg in , s bk f, and / t b a k f depending 
on whether both vertices are reconstructed with multiple tracks or not. The parameters for 
the background function Pbk g are determined using the AE-M^ sideband region for each 
decay mode. A MC study shows that the fraction of prompt component in the signal region 
is smaller (by ~10-50% depending on the decay mode) than that in the sideband region. 
We correct Pbkg for this effect. 

In the final fit to the events in the signal region, we determine simultaneously 12 param- 
eters: the B° and B~ lifetimes, 7 parameters for P rec and P asc , and 3 parameters for the 
outlier component. We find cr ol = 36^ ps, and / ol = (0.061°;^)% or (3.1 ± 0.4)% (multiple- 
or single-track case). The lifetime ratio, r = t b -/t b0 , is obtained by repeating the final fit 
after replacing r B - with rr-g . The fit yields r-g = 1.554 ±0.030 ps, t b - = 1.695 ±0.026 ps, 
and T B - /t b o = 1.091 ± 0.023. Figure @ shows the distributions of At for B° and B~ events 
in the signal region with the fitted curves superimposed. 

The systematic errors are summarized in Table |. The systematic error due to the IP 
constraint is estimated by varying (±10 /jm) the smearing used to account for the transverse 
B decay length. The IP profile is determined using two different methods and we find no 
change in the lifetime results. Possible systematic effects due to the track quality selection 
of the associated B decay vertices are studied by varying each criterion by 10%. The fit 
quality criterion for reconstructed vertices is varied from £ < 50 to £ < 200. We estimate 
the systematic uncertainty due to the maximum |At| requirement by varying the |At| range 
by ±30 ps and taking the maximum excursion to be the systematic error. We examine the 
uncertainty in the scale of At arising from the measurement error of the SVD sizes and 
thermal expansion during the operation, and find that its contribution to the lifetimes is 
negligibly small. AP-Mb c signal regions are varied by ±10 MeV for AE and ±3 MeV/c 2 
for M bc . The parameters determining / sig are varied by ±lcr to estimate the associated 
systematic error. The systematic error due to the modeling of P S j g is estimated by comparing 
the results with different P S j g parameterizations. The lifetime fit is repeated after varying the 
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FIG. 2: The At distributions of neutral (top) and charged (bottom) B meson pairs, with fitted 
curves. The dashed lines represent the sum of the background and outlier component, and the 
dotted lines represent the outlier component. 



TABLE I: Summary of systematic errors for neutral and charged B lifetimes, and their ratio. The 
errors are combined in quadrature. 





(P s ) 


Tb- (P s ) 


T B-/ T W 


IP constraint 


0.004 


0.003 


0.001 


Track selection 


0.006 


0.004 


0.001 


Vertex selection 


0.003 


0.002 


0.002 


At range 


0.003 


0.002 


0.001 


AE'-Mbc signal region 


0.003 


0.004 


0.003 


Signal fraction 


0.001 


0.001 


0.001 


R s i g parameterization 


0.008 


0.008 


Cancels 


i? np parameters 


0.006 


0.004 


0.006 


Background shape 


0.012 


0.007 


0.011 


MC statistics 


0.006 


0.007 


0.005 


Total 


0.019 


0.015 


0.014 
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R np parameters by ±2cr. The lifetime dependence on the B meson mass, which is the input 
for i?k, is measured by varying the mass by ±lcr from the world average value and found the 
differences to be negligible. The systematic error due to the background shape is estimated 
by varying its parameters by their errors. The possible bias in the fitting procedure and the 
effect of SVD alignment error are studied with MC samples; we find no bias and include the 
MC statistics systematic error. 

To conclude, we have presented new measurements of the B° and B~ meson lifetimes 
using 29.1 fb _1 of data collected with the Belle detector near the T(45) energy. Unbinned 
maximum likelihood fits to the distributions of the proper-time difference between the two 
B meson decays yield: 



These are currently the most precise measurements. A value of unity for tb-/t-^ is ruled 
out at a level greater than 3a and the measured value is consistent with the theoretical 
prediction [[J. 
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